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Abstract 

We study the (crcr), (ce), (ee) correlators in the 2d Ising model 
perturbed by a magnetic field. We compare the results of a set of 
high precision Montecarlo simulations with the predictions of two dif- 
ferent approximations: the Form Factor approach, based on the exact 
S-matrix description of the model, and a short distance perturbative 
expansion around the conformal point. Both methods give very good 
results, the first one performs better for distances larger than the 
correlation length, while the second one is more precise for distances 
smaller than the correlation length. In order to improve this agree- 
ment we extend the perturbative analysis to the second order in the 
derivatives of the OPE constants. 
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1 Introduction 



In these last years much progress has been done in the study of two di- 
mensional statistical systems in the neighbourhood of critical points. In the 
framework of quantum field theory these systems can be seen as Conformal 
Field Theories (CFTs) perturbed by some relevant operator. Since the sem- 
inal work of Belavin, Polyakov and Zamolodchikov |]J we have an almost 
complete understanding of CFT's (at least for the minimal models): we have 
complete lists of all the operators of the theories and explicit expressions for 
the correlators. However much less is known on their relevant perturbations. 
In some cases it has been possible to show that these perturbations give rise 
to integrable models @, |3|]. In these cases again we have a rather precise 
description of the theory. In particular it is possible to obtain the exact 
asymptotic expression for the large distance behaviour of the correlators ||]. 
From this information several important results (and in particular all the 
universal amplitude ratios) can be obtained. 

However in comparing with numerical simulations or with experiments 
one is often interested in the short distance behaviour of the correlators 
(short here means for distances smaller or equal than the correlation length) 
and this is not easily accessible in integrable systems. Moreover integrable 
perturbations represent only a small subset of the possible theories. For 
instance in the case of the Ising model both the purely thermal and the 
purely magnetic perturbations are integrable, but for any combination of 
them the exact integrability is lost. 

For these reasons, besides the S-matrix results, it is important to develop 
a perturbative approach well defined in the short distance regime of the the- 
ory and such that it does not rely on the exact integrability of the model. 
This is however a rather difficult task. In fact any naive perturbative expan- 
sion of the (massless) CFT along a relevant direction, is affected by infrared 
divergences (IR) and some non-trivial strategy is needed. 

Recently, in P, 0], a new approach has been proposed to overcome this 
difficulty (see |7||-|[L5|| for relevant related works and preexisting ideas). The 
method is based on Wilson's operator product expansion (OPE). Roughly 
speaking the main idea of this new approach is that the Wilson's coefficients 
of the OPE, being well defined at short distance, can be assumed to have a 
regular, IR safe, perturbative expansion with respect to the coupling. For 
this reason we shall refer to it in the following as the IRS (infrared safe) 
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perturbative approach. 

The main requirement of this IRS approach is the the knowledge of the 
Wilson coefficients (and their derivatives with respect to the perturbing cou- 
pling). For this reason it is particularly efficient if applied to perturbations 
of exactly solved theories like 2d critical CFT's (but the framework is quite 
general and in principle could be extended also to higher dimensions). 

The price one has to pay to control in this way the IR divergences is 
that one needs, as an external input information, the expectation values 
of the operators involved in the expansion. There are at this point two 
possibilities. The first one is to concentrate only on observables in which 
these expectation values exactly cancel. This is a small but very interesting 
subset of the informations that we can obtain with the IRS perturbation. 

The second possibility is to obtain the desired expectation values with 
some other method or extract them from numerical simulations (an inter- 
esting numerical approach to obtain these VEV is based on the Truncated 



Conformal Space technique, see |16], [17fl and references therein). 

From this point of view, the IRS approach becomes particularly pow- 
erful if applied to integrable perturbations, since in this case some of the 
expectation values can be deduced from the S-matrix of the model. 

The last step one has to face in comparing the results of the IRS method 
with simulations or experiments is the presence of a nonuniversal normaliza- 
tion factor between the operators in the continuum quantum field theory and 
their lattice discretizations. These normalizations (and the related normal- 
ization of the coupling of the perturbation) can be fixed if an exact solution 
of the lattice model exists at the critical point. Actually much less is needed. 
One only needs the exact expression (or even only its large distance asymp- 
totic form) of a correlator involving the operators in which we are interested. 
This makes the Ising model perturbed by a magnetic field a perfect candi- 
date for testing the IRS method. In fact it is well known that this model 
is exactly integrable ^ and all the amplitude ratios and expectation val- 
ues of the primary fields are known. Moreover the Ising model is exactly 
solvable at the critical point and the exact expression is known for several 
correlators [|18|, [L9| . 



In fact the IRS approach was successfully tested with the magnetic per- 
turbation of the Ising model in . The aim of this paper is to make a further 
step in this direction. In particular we have three goals: 
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a] Compare the results of the method with new high precision Montecarlo 
simulations so as to test the range of applicability of the method. 



b] Compare the IRS method with the results obtained in the S-matrix frame- 

work with the so called "form factor" (FF) approach. 

c] Show that it is possible to extend the analysis of f| |J to higher orders in 

the perturbative coupling and discuss the technical problems that one 
has to afford following this route. 

In particular we shall study in this paper, as an example, the second order 
term in the perturbative expansion of the (ee) correlator. The reason for this 
choice is that this correlator has a very peculiar behaviour since its first order 
correction turns out to be exactly zero, thus our second order calculation is 
mandatory if one is interested to study the influence of the magnetic field on 
the simple critical behaviour of the correlator. 

This paper is organized as follows. Sect. 2 is devoted to a general de- 
scription of the Ising model in a magnetic field both on the lattice and in the 
continuum. The aim of this section is to fix conventions and normalizations 
which will be useful in the following. In sect. 3 we shall briefly describe the 
IRS method, while in sect. 4 we shall extend it to second order derivatives of 
the magnetic field. In sect. 5 we shall briefly describe our Montecarlo simu- 
lation while in sect. 6 we shall compare the results of our simulations with 
IRS and FF predictions. Finally sect. 7 will be devoted to some concluding 
remarks. The details of the calculation of the second order derivative of the 
Wilson coefficient are collected in the Appendix. We have reported in three 
tables at the end of the paper a sample of the results of our simulations. 

2 Ising model in a magnetic field. 

The continuum theory, which is the starting point of the IRS expansion is 
given by the action: 



where Ao is the action of the conformal field theory which describes the 
Ising model at the critical point. Let us start our analysis by looking in detail 
at this CFT. 




(1) 
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2.1 The Ising model at the critical point 

The Ising model at the critical point is described by the unitary minimal 
CFT with central charge c = 1/2. It contains three conformal families whose 
primary fields 1, a, e have scaling dimensions x — 0, 1/8, 1 respectively The 
fusion rule algebra is 

[#] = [1] 

[a][e] = [a] (2) 
MM = [1] + [e]. 

Once the operator content is known, the only remaining information 
which is needed to completely identify the theory are the OPE constants. 
The OPE algebra is defined as 

^W^(o) = E4 fc} W<f w (o) (3) 

where with the notation {k} we mean that the sum runs over all the fields 
of the conformal family [k]. The structure functions CfAr) are c-number 
functions of r which must be single valued in order to take into account 
locality. In the large r limit they decay with a power like behaviour 

C*(r) ~ \r\~ dim{c ^ (4) 

whose amplitude is given by 

C£ = lim C£(r) \r\ dim ^\ (5) 

The actual value of these constants depends on the normalization of the 
fields, which can be chosen by fixing the long distance behaviour of, for 
instance, the aa and ee correlators. In this paper we follow the commonly 
adopted convention which isQ: 

(a(x)cr(O)) = r , \x\ — > oo (6) 

\x\ 4 



<e(x)e(0)) = -j— j2 , \x\^ oo. (7) 
\xr 



1 Notice the change of normalization with respect to ||. The e operator of the present 
paper corresponds to 2it times that of ref. || . 
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With these conventions we have, for the structure constants among primary 
fields 

= <% = q CT = o (8) 

= = Cl = = 1 (9) 

and 

= c^ a = \. (io) 

2.2 The Ising model in a magnetic field 

If we switch on the magnetic field h, the structure functions acquire a h 
dependence so that we have in general 

*i(r)*i(0) = £4* } (M*{*}(0) ■ (11) 

{fe} 

Also the mean values of the a and e operators acquire a dependence on 
h. Standard renormalization group arguments allow one to relate this h 
dependence to the scaling dimensions of the operators of the theory and lead 
to the following expressions: 

(a} h = A a hTs+... (12) 
(e) h = AM + - (13) 



The exact value of the two constants A a and A e can be found in []20[ and pT| 
respectively 

2C 2 

A a = — , 5 = 1.27758227.. , (14) 

15 (sin f + sin f + sin § ' v ; 



with 



and 



r (I) r i l r (J) r> (i) 1 



A e = 2.00314... . (16) 

Notice however that these amplitudes are not universal. They depend on the 
details of the regularization scheme. Thus some further work is needed to 
obtain their value on the lattice. 
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2.3 The lattice model 

The lattice version of the above model is defined by the following partition 
function: 

Z= /^W)^^^" 1 (17) 

(Ti=±l 

where the notation denotes nearest neighbour sites in the lattice which 
we assume to be a two dimensional square lattice of size L. In order to select 
only the magnetic perturbation, (3 must be fixed to its critical value: 

(3 = p c = \ log {V2 + 1) = 0.4406868... 
finally by defining hi = f3 c H we find 

Z = e pc ^^) aiaj+hl ^ iai (18) 

<Ti=±l 

In the following we shall denote the lattice discretization of the operators er, 
e with the index /. The magnetization M(h) is defined as usual: 

M(h) = ~Qog Z)\ p=Pc = (^E^)- (19) 

l i 

where N = L 2 denotes the number of sites of the lattice. This result suggests 
the following definition for the lattice discretization of o 

i 

so that the mean value of ai coincides with M(h): 

(ai) = M(h) (21) 
Similarly, we define the internal energy as: 



E M=^A^ z ^=^r' } (22) 
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For the energy operator one must also take into account the presence of an 
additional bulk contribution at the critical point. This constant can be easily 
evaluated (for instance by using Kramer- Wannier duality) to be e = \- 
This result suggests, for the lattice discretization of e, the following definition 

£ < S 2lv|> ff '-72 (23) 
so that the mean value of e; coincides with the singular part of E(h): 

(ei) = E(h) - -L (24) 
According to the above discussion we expect: 

{n) h = A l a hp + ... (25) 

(e l ) h = A[hp + ... (26) 

where the lattice amplitudes A l a , A\ are different from the corresponding 
amplitudes evaluated in the continuum. 

In order to relate the lattice results with the continuum ones, we must 
fix the relative normalizations of a versus o\ , e versus e\ and h versus 
The simplest way to do this is to look at the analogous of eq.(H,|7|) a t the 
critical point (namely for hi = 0) [J22|J . From the exact solution of the Ising 
model |L8| we know that 

Mh=o = jj^j-, (27) 



where denotes the distance on the lattice between the sites i and j and 



we know from |19[ that: 

Rl = e 3 «'(- 1 )2 5 / 24 = 0.70338... (28) 
By comparing this result with eq.(^]) we find 

0l = R aa = 0.83868... a (29) 

2 This essentially amounts to measure all the quantities in units of the lattice spacing. 
For this reason we can fix in the following the lattice spacing to 1 and neglect it. 
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From this we can also obtain the normalization of the lattice magnetic 
field which must exactly compensate that of the spin operator in the pertur- 
bation term ha. We find: 



hi = {R^h = 1.1923.. .h (30) 

Combining these two results we obtain the value in lattice units of the 
constant A a 

A l a = (R a ) mllb A a = 1.058... (31) 

In the case of the energy operator the connected correlator on the lattice, at 
hi = and for any value of /3, has the following expression [Q: 

2 



te(0)ei(r))« 



71 



Kl(5r)-Kl(5r) 



(32) 



where K and K\ are modified Bessel functions, 5 is a parameter related to 
the reduced temperature, defined as 



5 = A\(3-(3 C 



(33) 



and with the index c we denote the connected correlator (notice that thanks 
to the definition (p3|) no disconnected part must be subtracted at the critical 
point and the index c becomes redundant). This expression has a finite value 
in the 5 — ► limit (namely at the critical point). In fact the Bessel functions 
difference can be expanded in the small argument limit as 



Kl(5r) - Kl{8r 
thus giving, exactly at the critical point: 



(5r) 



+ 



(jrr)' 



Hence R e = By comparing this result with eq. (0) we find 

e 



R e e = - 

7T 



and from this we obtain the expression in lattice units of A £ 

A\ = {R a ) 8 / 15 {R e )A £ = 0.58051... 



(34) 



(35) 



(36) 



(37) 
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2.4 Correlators 



In the remaining part of this paper we shall be mainly interested in the 
dependence on the external magnetic field of the following correlators: 

G(t,c = (<K0)<r(r)) (38) 
G t , = (e(O)e(r)) (39) 
G a , e = k(a{0)e{r)) (40) 

where k = sign(h). We already know the behaviour at the critical point of 
the first two of them, which is given by eq. (0) in the continuum (or 

equivalently eq.(^), ( pop on the lattice), while the OPE constants reported 
in eq.(||) immediately tell us that (e(O)cx(r)) = 0. 

For small values of h we may expect to add to these results correction 
terms functions of h and r. However standard renormalization group argu- 
ments show that these two variables are actually related and that there is a 
natural scaling variable which describes the short distance expansion of these 
correlators in a magnetic field which is t = \h\ |r| 15 / 8 . In order to obtain an 
explicit expansion in powers of t we must absorb the scaling dimensions of 
the various operators in the expansion []. To this end let us define 

= {a{<S)a{r))\r\ x l A (41) 

F^= (e(0)e(r))|r| 2 (42) 
F CTi£ EEA;(a(0) e (r))|r| 9 / 8 (43) 

where k = sign(h). 

The powers which appear in the t expansion of the functions F can be 
immediately deduced from the analysis of the OPE via the IRS method, that 
will be described in the following section. 

3 On the contrary, in the large distance regime where one may use the predictions 
obtained with the form factor approach, the natural normalization is that of the 
functions defined above. 
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3 The infrared safe approach 



The goal of the method presented in Ref. [|5| is to obtain informations about 
the short distance behavior of a conformal field theory perturbed by relevant 
operators. 

The general idea behind this approach, (for preexisting ideas see || 
[X0| , PI , |T2| , |13] ) is the fact that Wilson Coefficients, being short distance 
objects, can be taken to have a regular, IR safe, perturbative expansion with 
respect to the coupling. This OPE approach leaves unfixed some constants 
that parameterize the vacuum expectation values of operators that appear 
in conformal field theory. 

In || it was found that the correlators of the perturbed CFT are given 
in terms of the derivatives of the Wilson coefficients (calculated at h — 
point). To be precise, they appear in the following way 



(Mr)^0))h = ^c!f } (h f r)^ w {0)) h 
{k} 



E 

{k} 



4* } (0,r)+^C#>(0,r) fc+Ia£C#>(0,r) h 2 + 



<*{*}(0))*. 



It was also shown that a general formula could be written for the n — th 
derivative of the Wilson Coefficients with respect to h. Here we will write 
only the first and the second order derivatives for the Wilson Coefficients, 



and 



fd 2 z([a z U a <P b -jZC b aia2 <P b \x R }) (44) 



= J'd 2 z fd 2 z'([a z a z , ^ b -j2C b aia2 <j )b S] jX R }} + 
+ it,d h C b aia2 fd 2 z([a z( f> b X R }). (45) 
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The general structure of this formula is a sum of the "naive" perturbative 
term plus n (for the n — th order coefficient) infrared counterterms. The 
asterisk reminds that the sum on the counterterms is truncated and is per- 
formed up to a given infrared dimension (see again ||). 

The construction of the IRS expansion requires two steps. 

• First, one must select by using the OPE rules which operators can ap- 
pear in the various expansions, identify their scaling dimensions, select 
the dominant ones and give their expression in terms of the structure 
constants and of their derivatives. 

• Second, one must evaluate the derivatives of the structure constants by 
using eq.fl44|) or ([45) to reduce them to suitable integrals over correla- 



tors evaluated at the critical point. This allows in principle to complete 
the analysis, since the explicit form of all possible critical correlators is 
known. However in general these integrals are highly non trivial and 
their evaluation represents the major problem of the whole approach. 

These two steps where performed in |J for all the terms in the expansion 
involving at most first order derivatives of the structure constants. This 
allows to obtain the first three terms in the expansion of the (ere) and (era) 
correlators (which are reported for completeness at the end of this section). 
On the contrary for the (ee) in this way one can only obtain the first two 
terms of the expansion. Moreover one can verify by an explicit calculation 
that the second of them is identically zero (||). Thus in the (ee) correlator, 
in order to reach the first non trivial correction to scaling, it is mandatory 
to extend the analysis of || and to deal with second order derivatives of 
the Wilson coefficients. We shall address this problem in the next section. 
In particular, in sect. 4.1 we shall discuss the first step of the IRS analysis, 
and select among the possible candidates the one with the lowest power of 
t which, as anticipated, turns out to involve a second order derivative of a 
structure constant. Then in sect. 4.2 (and in the Appendix) we shall explicitly 
evaluate this contribution. Let us conclude this section by listing for all the 
three correlators the first three terms of the IRS expansion 

F a ,a = Bl + B*„1*fl* + Blt 1 ^ + 0(t 2 ) (46) 
F e , € = B l te + Bit™' 15 + Bit 2 + 0(t 32 / 15 ) (47) 
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kF^ = Bit 1 ' 15 + B 2 ae t + B 3 a€ t 23 / 15 + 0(t 31 / 15 ) (48) 



where the coefficients B^ are 



Bl 




ci 


Bl 






Bl 






Bl 




ci 


Bl 




A a d h C- 


Bl 




\dtCl 


Bl 




A*Cj e 


Bl 




d h Cl t 


Bl 






given 


by 


(g4j) and 



the derivatives of C% h ^. are given by (|4]) and the notation d^C^. is the 
extension of the definition given in eq.@ to the derivatives of the Wilson 
coefficients. The first order derivatives have been calculated in || we report 
here for completeness their numerical value. Notice a slight change with 
respect to || due to the different choice of normalizations for the e operator 
(the e of the present paper corresponds to 2ir times that of ref. ||) 



d h C% a = -0.40374 

d h Cg = 

dlCl e = 3.29627 

d h C% e = -0.90900 

The second order derivative which appears in last one in the (ee) correlator 
requires a more involved calculation which we shall discuss in the next section 
and in the Appendix. 



4 Second order corrections 
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4.1 Dimensional analysis 

To estimate the higher order corrections to (ee) we must analyse two kinds 
of possible contributions. 

• The expectation values of secondary operators multiplied by the Wilson 
coefficients and their first derivatives. 

• The second derivatives of Wilson coefficients. 

We would like to understand which is the most important of these terms. 
In the Ising model there are two secondary operators at first level, obtained 
by acting on e and a with the Virasoro generator L_i and its hermitian 
conjugate (the action of L_i on 1 gives 0). We start by considering 

e 1 = L_iZ_ ie (49) 

and 

a 1 = L_xL-\(J (50) 

where are Virasoro generators. It is clear that the expectation value 

of this kind of operators is zero being total derivatives. So let us go to second 
level of the algebra. There are two possible terms: 

L 2 1 Z 2 x 0, L_ 2 I_ 2 (51) 

where is a generic primary field. In this situation we have to consider also 
the identity operator. 

In the identity sector the only contribution is given by 

Tf = L„ 2 L„ 2 1 (52) 

i.e. the energy-momentum tensor. Again a simple dimensional analysis shows 
that 

dimTf = 4, (TT) h = A TT \h\ 32/15 (53) 

giving 

A TT Gg t 32/15 . (54) 

It is clear that the terms containing secondary operators (of second level) of 
a and e are of higher order in t and will not be considered here. 
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A second possible contribution is given by the higher order derivative 

b 

= f d 2 z f d 2 z'([a z a z , I^M - £ C h aia2 ^ X R \) + 

+ E^a 1B2 fd 2 z([a z <f> b X R }). (55) 
b J 

Let us fix X R = 1. An elementary computation shows that the series is 
truncated and only those operators having x b < ^ appear in it. 
It follows that 



h " = 

= fd 2 z f ' d 2 z'(o z o z ,{e r e Q - C] € )) + 
+ d h C: e fd 2 z(a z a ) (56) 
but dhC° t = 0, and we can say that 

dlCl = f d 2 z J' d 2 z'(a z a z/ (e r e - C} e )). (57) 

Again from dimensional analysis we get 

. _ 14 

dim 9^ = - (58) 
and the contribution to (ee) of the second order derivative is given by 

\ dlCl t 2 . (59) 

It is also clear that, being X R = 1 the lowest dimension operator, deriva- 
tives of Wilson coefficients relative to a and e will give terms with a higher 
power in t. 

Let us write finally the perturbative expansion of this correlator 

Fee = Cf e + \d 2 & e t 2 + 0{e 2 ^). (60) 
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4.2 The Wilson derivative 



Let us remember that 



d 2 h Ci = f d 2 z f d 2 z'(a z a z ,(e r e - (61) 

where (cr(zi)a(z2)e(z 3 )e(z4)) denotes the correlator at the critical point which 
can be written as 



(^w*)«(*w«o> = :" raT "r: ff ° ' (62) 



kl2( 2 32 + z 42, 


— 2^32^42 


2 


4 1 Z42Z32Z41Z31 \ 


^43 


2 


^12 


1/4 



By fixing the values of z\ = z, Z2 = w, Z4 = and by rescaling r we can 
choose z 3 = 1, we get 

^1 _ f d 2 r d 2 z \Z(1-W)+W(1-Z)\ 2 

where the dots indicate the counterterms. 

The explicit calculation of this integral can be done using a technique 



developed by Mathur, [23]. The general idea behind this approach is to 
factorize the integral in a holomorphic and antiholomorphic part using Stokes 
theorem. The calculation is reported in the Appendix. After this calculation, 
in order to get rid of the infrared cutoff we perform a Mellin transform of the 
integral and the infrared counterterm (see for more details). In this we 
we end up with a finite result when the infrared cutoff goes to infinity. The 
final result is 

did = 97.5936.... (64) 

Let us stress that the techniques that we have discussed in this section 
can be extended to any order in the derivatives of the Wilson coefficients. 
This is a rather important observation, since it allows, in principle, to study 
the IRS corrections, in a consistent way, to any given order in t. 



5 The Montecarlo simulation 



It has been recently shown [25] that in the case of the 2d Ising model in a 
magnetic field, algorithms based on the exact (or approximate) diagonaliza- 
tion of the transfer matrix are much more effective than standard Montecarlo 
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simulations. In particular this is true for all possible observables involved in 
the large distance behaviour of the model. The only exception is represented 
by the short distance behaviour of the point-point correlators which is the 
subject of the present paper. In fact in order to reach lattices as large as 
possible in the transfer-matrix programs discussed in [^5] only zero momen- 
tum projected observables could be studied, while we are instead interested 
in point-point correlators. Moreover, we need to have a window as large as 
possible between the region (few lattice spacings) dominated by the lattice 
artifacts and the correlation length. This windows shrinks to zero in the 
transfer matrix approach where only small values of the correlation length 
can be studied. 

For these reasons we decided to perform our tests with standard Monte- 
carlo simulations. We used a Swendsen-Wang type algorithm, modified so 
as to take into account the presence of an external magnetic field. For a 
detailed description of the algorithm see for instance [21| ^7 . 



5.1 Finite size effects. 

As a preliminary test we performed a simulation at hi = 4.4069 x 10~ 4 (which 
corresponds to H = 0.001) with lattice size L = 128 which exactly coincides 



with one of the simulations reported in |26] and found results in complete 
agreement with those quoted in [p6fl . Then for the same value of hi we 
performed a set of high precision simulations varying the lattice size so as to 
check the presence of possible finite size effects. In particular we compared 
our estimates of the mean magnetization, susceptibility and internal energy 
with the known exact results, extrapolated at the value of hi at which we 
performed the simulations 0. 

The comparison is reported in tab.l. It turns out that lattice sizes at 



4 If one is interested in a high precision comparison, also the contribution of secondary 
fields should be taken into account in extracting these exact estimates. The amplitude of 
some of these secondary fields have been evaluated numerically in |^5| . In the case of M 
and Xj f° r the values of hi in which we are interested, the contributions of the secondary 
fields are strictly smaller than the statistical errors of the results of our simulations and 
hence can be neglected in the comparison. On the contrary for the internal energy it turns 
out that the amplitude of the first correction is rather large (see [^5| for details) and must 
be taken into account. In fact, if we would neglect it, instead of the value reported in 
tab.l we would find E = 0.71652 in clear disagreement with the Montecarlo results. This 
represents a non trivial test of the results of |23] . 
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Table 1: Finite size effects at hi = 4.4069 x 10~ 4 . tn the hrst column we report 
the lattice sizes used in the simulations, in the remaining three columns the 
mean values of the magnetization, susceptibility and internal energy. In the 
last row we report the exact results obtained by using the known values of 
the amplitudes for these quantities. 



L 


M 


X 


E 


120 


0.63110(16) 


113(1) 


0.71638(3) 


140 


0.63247(16) 


98.4(6) 


0.71645(3) 


160 


0.63245(14) 


96.4(4) 


0.71639(3) 


200 


0.63255(11) 


95.4(3) 


0.71643(3) 


exact 


0.63260 


95.7 


0.71642 



least larger than 12 times the correlation length are needed to be sure that 
finite size effects are under control (with this we mean that the systematic 
errors induced by the finite size of the lattice are smaller than the statistical 
errors of the simulations and can be neglected). A side consequence of this 
observation is that the simulations reported in |26j are indeed affected by 
rather large finite size effects. 

It is interesting to notice that the magnetic observables are more affected 
by finite size effects than the thermal one. As one can easily expect the 
largest corrections appear in the case of the susceptibility. 

5.2 The simulation 

Once we were sure to have finite size effects under control we performed a 
set of high precision simulations of the model for three different values of the 
magnetic field. 

An important quantity to understand the range of validity of the IRS 
approximation is the correlation length. Roughly speaking we expect that 
the IRS results should give a reasonable approximation for distances equal or 
smaller than the correlation length, while above it the form factor approach 
should give results of better quality. For this reason it is important to have 
a good estimate of £. This can be easily obtained from the knowledge of the 
spectrum of the theory. We find, in lattice units: 



17 



£(hi) = 0.24935. ..fy 15 (65) 

(see p5| for details on the continuum to lattice conversion of £). In tab. 
2 we have reported the expected values of £ in our cases. 



Table 2: Values of the correlations lengths for the three choices of hi. 



hi 




4.4069 x 10" 4 
2.2034 x 10" 4 
1.1017 x 10" 4 


15.4 
22.4 
32.2 



For all the values of h\ that we simulated, we studied the three cor- 
relators: (cr(O)cr(r)), (er(O)e(r)) and (e(O)e(r)), for r = 1, ■ ■ ■ , L max , where 
the maximum distance L max was chosen to be roughly twice the correlation 
length. In this way we can test our results also in the large distance regime, 
where predictions from the form factor approach are expected to give very 
precise estimates for the correlators. Notice that when studying the large dis- 
tance behaviour of correlators one is usually interested in the zero momentum 
projection of the connected part of the correlator. On the contrary in the 
present case we are interested in the point-point correlators without mean 
value subtraction or zero momentum projections. This must be taken into 
account when comparing the data with those obtained with the form factor 
approach. Some informations on the simulations are reported in tab. 3. 



Table 3: Some informations on the simulations. L max denotes the maximum 
distance at which the correlators have been evaluated, its value almost co- 
incides with twice the correlation length. L denotes the lattice size, hi the 
magnetic held. In the third column we have reported the number of measures 
while in the fourth column we have reported the number of SW sweeps which 
separates two measures. 



hi 


L 


measures 


sweep / measures 


Lmax 


4.4069 x 10~ 4 


200 


4 x 10 5 


5 


30 


2.2034 x 10~ 4 


300 


2 x 10 5 


5 


45 


1.1017 x 10~ 4 


400 


1 x 10 5 


10 


65 
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We report an example of our results (for the value hi = 4.4069 x 1CT 4 ) in 
the first columns of tabs. 7, 8 and 9. The quoted errors have been obtained 
with a standard jacknife method. 



6 Discussion of the results 

In fig.s 1-8 and tab. 7, 8 and 9 we compare our estimates for the correlators 
with the IRS and form factor predictions. For completeness we briefly recall 
here the form factor results (see ^9[ for details) and give the numerical 
values (once all the conversion factors are taken into account) of the constants 
in the IRS approach. 



6.1 Form factor results 

The scattering theory which describes the scaling limit of the Ising Model in 
a magnetic field |§ contains eight different species of self-conjugated particles 
A a , a = 1, ... ,8 with masses 

m 2 = 2m x cos - = (1.6180339887..) m x , 
5 

m 3 = 2m 1 cos— = (1.9890437907.. )m u 
30 

7vr 

m 4 = 2m 2 cos— = (2.4048671724..) mi, 
30 

2-7T 

m 5 = 2m 2 cos— = (2.9562952015.. )m u (66) 

m 6 = 2m 2 cos— = (3. 2183404585.. )m 1 , 
30 

7T 7lT 

m7 = 4m 2 cos — cos — = (3.8911568233..) mi , 
5 30 V ' ' 

m 8 = 4m 2 cos — cos — = (4.7833861168..) mi 
5 15 

mi(hi) denotes the overall mass scale and coincides with l/£, hence its value 
in lattice units is 

mx{hi) = 4.0104... hf . (67) 

From the knowledge of the masses and of the form factors it is possible 
to obtain a large distance approximation for the correlators by constructing 



19 



a spectral sum over a complete set of intermediate states. We thus find for 
any pair of local operators $1 and $2: 



GWx) = ($i(z)$ 2 (0)) 



n=0 

X 



I 33 ! 52fe=i mfcCosh(9fc ^ ^gg-j 



where the form factors a {6\, . . . , 9„) are the matrix elements of the local 
operator on the asymptotic states A a , i.e they are defined as 



Fl..aS0u . . . A) = <0|$(0) 14^(01) • • • 4J0„)> • (69) 

The important point is that these form factors can be exactly computed 
in the integrable models once the S'-matrix is known. 

It is natural to organize the above expansion setting a reference value m r 
and keeping in the spectral sum only the states with a mass smaller than 
m r . Looking at eq.(p§D we see that we must expect, as a consequence of this 
truncation, a systematic error in the approximation of the order O (e~ mrX ). 
Up to m r = 2mx (where mi is the fundamental mass of the model) only 
single particle states survive in the sum and eq.fl68|) greatly simplifies. In 
particular looking at eq.(66) we see that only the first three states (which are 
the only ones below the pair production threshold) survive. Eq. ( |B8D becomes 
in this case, (choosing for instance the spin-spin correlator) 

GW(r) ~ M 2 (^)(l + E ^# (w)) (70) 

i=i * 

where K (x) is the modified Bessel function and Ff denotes the overlap 
(measured in units of the magnetization) of the i th state with the a operator. 
M(hi) denotes the magnetization and its hi dependence is given in ([L2|, |31|) . 
Similarly the spin-energy and energy-energy correlators are given by 

3 per* rps 

kG ae {r) ~ M{hi)E(hi)(l + £ — — *-#o(w)) (71) 

3 I I ^ 

G et (r) ~ E^h^l + £ ^i^oKr)) (72) 
i=i n 
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The constants F° and F[ have been evaluated in [B3, Their value is 
reported in tab. 4 and 5. 

One can systematically improve the approximation by setting higher val- 
ues of m r . For instance, for m r = 3 one must keep into account the first five 
single particle form factors together with the two-particle ones involving the 
(1,1), (1,2) and (1,3) pairs, and so on. By using the results of [|^, also 
these multiparticle form factors can be evaluated exactly. 



Table 4: Overlap amplitudes for the spin operator. 



iff = 


-0.64090211 


n = 


0.33867436 


n = 


-0.18662854 


n = 


0.14277176 


n = 


0.06032607 


F a — 


-0.04338937 


F a — 


0.01642569 


F a — 


-0.00303607 



Table 5: Overlap amplitudes for the energy operator. 



Ff = 


-3.70658437 


n = 


3.42228876 


n = 


-2.38433446 


n = 


2.26840624 


F§ = 


1.21338371 


Fl = 


-0.96176431 


F e — 


0.45230320 


F £ — 


-0.10584899 



As discussed above, in eq.s (f70[) , (|7T|) and (|72| ) we expect systematic errors 
of the order O (e~ 2mir ). For distances larger than the correlation length these 
deviations are very small but become increasingly relevant as the correlation 
length is approached. It would be important to have an estimate of the 
magnitude of these corrections. From this point of view the present case is a 
perfect laboratory since we know that all other possible sources of systematic 
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errors are under control. We report an example of the results obtained using 
eq.s (0), (|TT]) and fl7|) (for the value hi = 4.4069 x 10" 4 ) in the last column 
of tabs. 7, 8 and 9. In order to gain further insight on the convergence of 
the approximation we have also evaluated the contribution of the first eight 
terms of the spectral series (i.e. with m r = 3m>i). We report in fig. 10 the 
result of this analysis, together with those with m r = 1m,\ and m r = mi for 
comparison, in the particular case of the aa correlator. 

6.2 IRS approach 

By using the results of sect. 3 and 4 and the exact knowledge of the constants 
R a and R e one can easily write the constants in lattice units. They are 
reported in tab. 60. 

We report an example of our results (obtained by plugging the values of 
tab. 6 in eqs.(f46|), (f47|) and (gHD in the particular case hi = 4.4069 x 10" 4 ), in 
the second column of tabs. 7, 8 and 9 . 

Table 6: Coefficients of the IRS expansion in lattice units. 





= 0.7034... 


BL 


= 0.6414... 


BL 


= -0.3007... 


Bl 


= 0.1013... 


Bl 


= 


Bl 


= 3.4776... 


BL 


= 0.1685... 


Bl 


= 0.7380... 


BL 


= -0.3712... 



6.3 Comparison with MC results 



5 In doing this conversion one must also take into account the h factor contained in t. 
To stress this fact we introduce in analogy to hi a new scaling variable ti = \h[\ |r| 15 / 8 . In 
the coefficients reported in tab. 6 the conversion from t to ti has already been taken into 
account, hence they refer to the IRS expansion in powers of ti of the correlators 
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6.3.1 Lattice artifacts 

It is interesting to see that lattice artifacts are confined to a remarkably small 
region of few lattice spacings, which shows a negligible dependence on the 
magnetic field or on the type of correlator. Since the lattice artifacts decrease 
so quickly it is rather easy to find where does the region of applicability of 
the IRS results starts, by simply looking at the distance L min where for the 
first time the IRS prediction becomes compatible (within the errors) with the 
MC data or, if this never happens, looking at the location of the minimum 
difference between IRS predictions and the MC simulations Q It turns out 
that for all correlators and hi values L min ranges between 7 and 9 lattice 
spacings. This tells us that, at least in the case of the Ising model perturbed 
by a magnetic field, the IRS method has a large window of applicability, 
which becomes larger and larger as the critical point is approached. This is 
well exemplified by fig. 6 and 7 where the difference between MC data and 
IRS predictions is plotted, for the aa correlator in the short range region, 
first in units of the lattice spacing and then in units of the correlation length. 
It would be interesting to test if this behaviour also holds for other models 
or for different realizations of this one. 

6.3.2 hi dependence 

In the range L min < r < £ the agreement between IRS predictions and 
MC results is always very good. In particular, it seems that the method 
reaches its better results in the case of ea correlator. As expected, the IRS 
approximation becomes better and better as we approach the critical point 
(see fig.s 3, 4 and 5). First because the range of validity becomes larger 
and second because the systematic deviations due to the terms neglected in 
the expansion, which are proportional to higher powers of hi become less 
important. In particular for the ea correlator at hi = 1.1017 x 10~ 4 there 
is a wide range (more than 40 lattice spacings) in which the IRS prediction 
coincides with the MC results within the errors (see fig. 8). 



6 We may expect that higher orders in the t expansion improve the large distance 
behaviour of the IRS results, but they give a negligible contribution around L m i n where 
the discrepancy between Montecarlo data and IRS predictions is completely dominated by 
lattice artifacts. 
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6.3.3 IRS versus FF. 



Looking at tables 7, 8, 9 and at the fig.s 1-5 we see that, as expected, the 
FF approach performs better than the IRS one for distances larger than the 
correlation length and that the opposite is true for distances smaller than £. 
It is interesting to see that for distances of the order of the correlation length 
the IRS and FF methods give comparable performances. Some interesting 
informations on the systematic errors involved in the two approximations can 
be extracted from the data. 

1] While the systematic errors in the IRS approach have a polynomial be- 
haviour in the distance r (which is contained in ti), those of the FF 
approach have an exponential behaviour. This is clearly visible in fig. 2 
where the deviations are plotted as a function of r for the correlator 
G aa at hi = 4.4069 x 10~ 4 and in fig.3-5 where they are plotted for all 
the correlators and all the values of hi as functions of r/£. This makes 
the IRS method still reasonably reliable even for distances twice the 
correlation length. 

2] We may obtain a rough estimate of the magnitude of the systematic er- 
rors involved in the IRS approach with the following argument. Since 
ti = \hi\ |r| 15//8 , looking at eq.(^5|) we see that for distances of the 
order of the correlation length we have ti ~ 0.06. Depending on 
the correlator chosen, we would expect deviations of the order 0(tf), 
0(tf 1//15 ), 0(tf 2 ^ 15 ). Since in general the -E>| $ constants are of order 
unity, this amounts to an expected deviation for the F&& functions of 
order SF ~ 0.004. This expectation is in good agreement with the 
values of SF obtained by comparing the MC and IRS estimates at the 
distance r = £. (see for instance the data reported in tab.s 7,8 and 
9. In using these data one must take into account the normalization 
between and Gq,^ functions). 

3] A similar analysis can be performed in the case of the FF method. In 
this case we know that the systematic errors are of order 0(e~ mrr ). 
We shall further discuss them in sect. 6. 3. 5 below. In the large distance 
regime r > 1.5 £ the performances of the FF approach are very good. 
For instance, in this region, for the lowest value of h that we studied: 
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hi = 1.1017 x 10 4 the FF predictions for the ee correlator coincide 
with the MC results within the errors (see fig. 5). 

6.3.4 Convergence of the IRS expansion 

It is interesting to study the convergence properties of the IRS method, i.e. 
to see if the agreement with the Montecarlo data improves as higher terms 
are added in the expansion. 

• (eer) and {era). 

In these two cases the agreement improves as new terms are added in 
the expansion. This is clearly visible in fig. 8 where we have plotted 
the difference between the MC data for the correlator (ea) at hi = 
1.1017 x 10~ 4 and the IRS results with one (pluses), two (crosses) and 
three (diamonds) terms in the expansion. The analogous plot for the 
(a a) correlator shows exactly the same behaviour. 

• (ee). 

In this case we find exactly the opposite behaviour. As it is shown 
in fig. 9 the new coefficient that we evaluated does not improve the 
agreement with the Montecarlo data which is, by the way, impressively 
good already with the simple zero order contribution to the correlator. 
We see two possible reasons for this, rather unexpected, behaviour. 

1] As noticed in sect. 4, the next term in the perturbative expansion of 
the correlator has an exponent 32/15, which is very near to the one 
that we evaluated. We cannot evaluate this further contribution 
since it would require the knowledge of the expectation value of 
the TT operator. In principle this term could well compensate 
the deviation that we observe. 

2] This behaviour could be an indication of the bad convergence prop- 
erties of the IRS method. If this is the case it would be very 
interesting to test (in view of the good behaviour of the other two 
correlators) if this is a peculiar feature of the (ee) correlator or a 
feature of the expansion itself. This issue could be settled in prin- 
ciple by looking to the higher perturbative terms in the expansions 
of the two other correlators. We plan to address this point in a 
forthcoming paper. 
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6.3.5 Convergence of the FF expansion 

In order to study the convergence of the FF approximation we have compared 
the Montecarlo data in the case of the (aa) correlator at hi = 1.1017 x 10~ 4 
with the result of the FF approximation truncated at m r = m\, m r = 2mi, 
m r = 3mi respectively. This corresponds to take into account one, three and 
eight states respectively in the spectral sum. The results of the comparison 
are reported in fig. 10 . Looking at this figure one may see that as higher 
orders are added the approximation smoothly converges to the MC data. 
By comparing the three approximations one may get a perception of the 
convergence rate of the method. 



7 Concluding remarks 

In this paper we have compared the predictions of the IRS and FF approx- 
imations for the aa, ea and ee correlators with the results of a set of high 
precision MC simulations of the 2d Ising model perturbed by a magnetic field. 
To this end we have extended the IRS approach to second order derivatives 
of the structure constants. Our main results are: 

• Lattice artifacts are confined in a small region of few lattice spacings. 

• There is a wide region ranging from ~ 7 — 9 lattice spacings to the 
correlation length in which the MC data are in good agreement with 
the IRS results. 

• The agreement improves as the critical point is approached. 

• For distances smaller than £ the IRS gives a better approximation than 
the FF method, while the opposite is true for distances larger than £. 

The IRS method can be extended in principle to any order in the deriva- 



tives of the Wilson coefficients, by using the integration method of |24] and 
the technique of the Mellin transform. However in the case that we studied, 
i.e. the (ee) correlator, the IRS method turns out to show rather bad con- 
vergence properties. It remains an open problem to understand if this is a 
limit of the method itself or if it is a peculiar feature of the correlator that 
we have chosen. 
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It would be very interesting to extend this analysis to other models in 
this same universality class. In particular one could study the model recently 
introduced in |J{], [TTJ for which an exact bethe ansatz solution, out of the 
critical point exists. Another interesting application of the method would be 
the study of the correlators in the case of the most general perturbation of the 
Ising critical point (i.e. a mixed situation with both magnetic and thermal 
perturbations). In this case the exact integrability is lost but the IRS method 
is still valid and could give important informations on the behaviour of the 
correlators. In particular it would allow us to compare our approximation 



with the interesting results, directly obtained on the lattice, in |32 



Appendix 

The evaluation of the coefficient d\C\ t involves the calculation of the integral 

Z = J d 2 w\w\ e \l-w\ f w* r (l~w) s J d 2 z\z\ a \l-z\^z n {l- z)* m \z-w\' y (73) 

where n, m, r, s G N and a, j3, 7, e, / G R. 
It is useful to introduce the following theorem (see 
an integral of the form 

N 

I = / d 2 w £ f a (w)Q a pfp(w*) (74) 

where {f a (w)}^ =1N , {fp(w*)}p =liN are two sets of independent functions 
and Q a /3 is a constant matrix. Let us assume that fp(w*) e (fp(w))* have 
the same monodromies, in particular the two sets of functions f a (w) and 
g = (fp(w*))* must have the same branch points {wk}™=o such that 

= I wo I < \wx\ < ■ ■ ■ < \w m \ < |w m +i| = 00 (75) 

and they have to be analytic elsewhere. 

If we assume now that the matrix Q is invariant under the monodromy group 
action 

Q = MlQM*, Vfc (76) 



, [pi). Let us consider 
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where are the monodromy matrices of / and g related to the branch 
points Wk, it follows that we are able to express I in terms of one-dimensional 
integrals (see ]23j], |J for more details) 



k+lj 



1 - A4)- 1 ) Q 



k=l 



ad 



f (k) 



(77) 



where * is the transposition and 

t ik) _ 



dwf(w) 
dw*f(w*) 



(78) 



where Ck (Ck) are counter-clockwise (clockwise) circumferences enclosing all 
the branch points of modulus lower than Wk, starting at Wk + (infinitesimally 
over the cut at Wk) and ending at Wk_ (infinitesimally under the cut at Wk). 

Now we are able to evaluate both the z-plane and iy-plane integrations 
of ( ]T5D using the previous lemma. 
First, to perform the £-plane integration, we pose 



I z (w,w*) = I d 2 z\z\ a \l - zfz^l - z)* m \z - w\' f 

so the z-plane integration involves the following branch points 

z = 0, Zi = w, z 2 = 1, z 
Thus, the application of ( ff7|) gives 



oo. 



I z (w,w*) 



l r 



x l 1 12 x 2 



+ 



^ r 



T (2) T (2)^-(2) 
x l J 12 x 2 



where 



and 



T 



T (i) 

1 12 



(2) 
12 



1) 



e ivr(a+7)^ e i7r/3 _ 
^ e i7r(a+/3+7) _ X)( e i7r ( a +T) — 1) 



(79) 
(80) 

(81) 
(82) 
(83) 
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are the only non vanishing entries of the matrices 

T« = ((l-M 2 )- 1 -(l-Af 1 )- 1 )*g 
T® = ((l-M 00 )- 1 -(l-M 2 )- 1 ) < g. (84) 
This imply that we have to take in account only the integrals 

T (l) _ ,-wa -, x ,l+^/2+ 7 /2+n r (Q/ 2 + " + jWy/g + 1) 

11 ~ {e ~ 1)W r(a/2 + 7 /2 + 2 + n) ' 

■ F(-/3/2,a/2 + n + l;a/2 + 7/2 + 2+n;w); 

22 " (e r(a/2 + 7 /2 + 2) ' 

• F(-m-/3/2,a/2 + l;a/2 + 7/2 + 2;w*) (85) 

and 

T (2) ^-Mo+ZH-r) 1V VTO r(-q/2 - >g/2 - 7 /2 - n - l)r(/?/2 + 1) 
±l ~ {e 1){ } r(-a/2- 7 /2-n) 

■ F(-j/2, -a/2 - /3/2 - 7 /2 - n - 1; -a/2 - 7 /2 - n; w); 
t(2 ) (M a^) , , r(-a/2 - 0/2 - 7/2 - m - l)rQ3/2 + 1 + m) 
2 " ^ " } r(-a/2- 7 /2) 

• F(-7/2, -a/2 - /3/2 - 7/2 - m - 1; -a/2 - 7/2; w*). (86) 

Finally, putting all these relations in (|31~D, we can recover the wanted result 
for J 2 (w, w*). 

The w-plane integration is very similar to the previous one. Now we have 
to evaluate 

Z= J d 2 w\w\ e \l-w\ f w*'\l-w) s I z (w,w*) (87) 

which involves wq — 0, W\ — 1, = oo as branch points. Hence the 
solution is given by (ff7|), i.e. 



Z = - 

2 



(88) 



where 

rW = ((1 - Moo) -1 - (1 - Mo)- 1 )* Q. (89) 
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The contribution coming from is 

= (e me - 1) f 1 dz/i = (e me - 1) Ji 

JO 

= / e -i7r(e+a+ 7 ) _ ^ /" ^ = / -^(e+a+7) _ 1) J 2 
JO 

= (e" i7re - 1) / d^/i = (e" i7re - 1) Ji 



j(l) = ( e Me+«+7) _ !) Z* 1 ^*^ = ( e Me+«+7) _ !) J 2 (90) 

Jo 

that, in terms of generalized hypergeometric functions, becomes 

Ji = S(-a/2-/3/2- 7 /2-n-l,/3/2 + l)S(l + e/2,l + //2 + s) 

3^ 2 ( 2 71 ' ~ 2 ' 2' 2 ~1T~ + s )5 X ) 

J 2 = + 7 /2, 1 + a/2 + ra)B(2 + e/2 + a/2 + 7 /2 + n, 1 + //2 + s) 

^ , P ~ a + 7 ,„ „ a „ a + 7 

3 F 2 (-|,2 + -^ + n + e/2,l + - + n;2 + -^- + n,3 + 

a+7+e+f 
+ g + ^ + s); 1) 

J x = B(-a/2 - /3/2 - 7/2 - m - l,/3/2 + 1 + m)B(l + e/2 + /c, 1 + //2) 
3F 2 ( ~ a ~ 2 /3 ~ 7 -m-l,^,l + e/2 + fc;^ 2 ^,2 + ^ + fe);l) 

J 2 = fi(l + 7 /2,l + a/2)5(2 + e/2 + a/2 + 7 /2 + A;,l + //2) 

^,—13 a + 7 + e , a a + 7 
3 F 2 (^-m,2 + i + fc > i + _;2 + -^-I 3 + 



+ a + 7 : e + / +fr);i). 



(91) 



Thus the solution has the form 

.Z = t\\J\J\ + t\lJ\ J 2 + ^21^2^1 + ^22^2^2 (92) 

where the matrix elements are the following 

f n = A- 1 1 S(e/2) 1 S(/?/2) 1 S((a + /3 + 7 )/2)- 
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• [s{a/2)S{(3/2)S{f/2)S{{a + p + e + f + 2 7 )/2) + 
+ Sb/2)S((a + /3 + l)/2)S((a + e + f + j)/2)S((f + 7 + 0)/2)) 



t 22 = A- 1 S{a/2)S{ 1 /2)S{{a + e + 1 )/2)- 

's{*/2)S{p/2)S{{e + f)/2)S(l/2((3 + f + 7)) + 

+ S( 7 /2)5((a + /? + i)/2)S(f/2)S({e + f + 7 + 0)/2)) 

* 12 = t 12 = Ar 1 S(a/2)S(f3/2)S(e/2)S( , y/2) ■ 

■ 5((a + e + 7 )/2)S((a + /3 + 7 )/2)S((/3 + 7 )/2) 

with 

A = S((a + 7 )/2) 

\s(a/2)S(P/2)S((e + f)/2)S((a + (3 + e + f + 2 T )/2) + 

+ S(r//2)S((a + (3 + 7 )/2)5((a + e + / + 7 )/2)S((e + / + 7 + /?)/2) 
and S'(x) = sin(7rx). 

For all details on the calculation we refer to I03 . 
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Fig. 1: Comparison of the Montecarlo estimates (crosses), the IRS results (aster- 
isks) and the form factor results (pluses) for the correlator G aiJ at h\ = 4.4069 x 1CT 4 . 
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Fig. 2: Differences between MC estimates and IRS results (crosses) and between 
MC and form factor results (pluses) for the correlator G ea at hi = 1.1017 x 1CT 4 . 
In this figure, and in all the following ones the errors in the MC estimates are not 
reported since they are smaller than the symbol size. 
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Fig. 3: Differences between MC estimates and IRS results (crosses for hi = 
4.4069 x 1(T 4 , dotted squares for 2.2034 x 10" 4 and circles for 1.1017 x 10" 4 j and 
between MC and form factor results (pluses for hi = 4.4069 x 10~ 4 , diamonds for 
2.2034 x 10~ 4 and filled squares for 1.1017 x 10 4 Jdiamonds) for the correlator G acr . 
Distances are measured in units of the correlation length. 
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Fig. 4: Same as fig-3, but for the G ae correlator. Notice the different scale on the y 
axis 
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Fig. 5: Same as fig. 3, but for the G ee correlator. 
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Fig. 6: Differences between MC estimates and IRS results for hi = 4.4069 x 1CT 4 
(pluses), hi = 2.2034 x 10~ 4 (crosses) and hi = 1.1017 x 10 -4 (diamonds) for the ea 
correlator, in the short distance region dominated by the lattice artifacts. 
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Fig. 7: Same as in fig.6, but with distances measured in units of the correlation 
length. 
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Fig. 8: Difference between the MC data for the (eo) correlator at hi = 1.1017x 10~ 4 
and the IRS results with one (pluses), two (crosses) and three (diamonds) terms in 
the expansion. 
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Fig. 9: Difference between the MC data for the (ee) correlator at hi = 1.1017 x 10 4 
and the IRS results with one (pluses) and two (crosses) terms in the expansion. 
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Fig. 10: Difference between the MC data for the (era) correlator at hi = 1.1017 x 
10~ 4 and the FF results with one (pluses), three (crosses) and eight (diamonds) 
terms in the spectral series. 
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Table 7: Comparison of the Montecarlo estimates (second column), the IRS 
results (third column) and the form factor (FF) results (fourth column) for 
the correlator G aa at hi = 4.4069 x 1CT 4 . In the first column is reported the 
distance in lattice units. 
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Table 8: Comparison of the Montecarlo estimates (second column), the IRS 
results (third column) and the form factor (FF) results (fourth column) for 
the correlator G ee at hi = 4.4069 x 1CT 4 . In the first column is reported the 
distance in lattice units. 
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Table 9: Comparison of the Montecarlo estimates (second column), the IRS 
results (third column) and the form factor (FF) results (fourth column) for 
the correlator G ae at hi = 4.4069 x 1CT 4 . In the first column is reported the 
distance in lattice units. 
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